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ABSTRACT 

We have measured gas and stellar velocity dispersions in 5 circumnuclear starform- 
ing regions (CNSFRs) and the nucleus of the barred spiral galaxy NGC 3351. The 
stellar dispersions have been obtained from high resolution spectra of the CaT lines 
at AA8494, 8542, 8662 A, while the gas velocity dispersions have been measured by 
Gaussian fits to the H/3A4861 A line on high dispersion spectra. 

The CNSFRs, with sizes of about 100 to 150 pc in diameter, are seen to be com- 
posed of several individual star clusters with sizes between 1.7 and 4.9 pc on an HST 
image. Using the stellar velocity dispersions, we have derived dynamical masses for the 
entire starforming complexes and for the individual star clusters. Values of the stellar 
velocity dispersions are between 39 and 67kms^^. Dynamical masses for the whole 
CNSFRs are between 4.9 x 10^ and 4.3 x 10^ M© and between 1.8 and 8.7 x 10^ M© 
for the individual star clusters. 

Stellar and gas velocity dispersions are found to differ by about 20kms~^ with 
the H/3 lines being narrower than both the stellar lines and the [Oiii] A 5007 A lines. 
We have found indications for the presence of two different kinematical components in 
the ionised gas of the regions. The radial velocity curve shows deviation from circular 
motions for the ionised hydrogen consistent with its infall towards the central regions 
of the galaxy at a velocity of about 25kms~^. To disentangle the origin of these two 
components it will be necessary to map these regions with high spectral and spatial 
resolution and much better S/N in particular for the O^^ lines. 

Key words: galaxies: individual (NGC 3351) - galaxies: starburst - galaxies: kine- 
matics and dynamics galaxies: star clusters (ISM:) HII regions 



1 INTRODUCTION 

The bulges of some nearby spiral galaxies show intense star- 
forming regions located in a roughly annular pattern around 
their nuclei. At optical wavelengths, these circumnuclear 
starforming regions (CNSFRs) are easily observable rings. 
In the ultraviolet (UV), massive stars dominate the observed 
circumnuclear emission even in the presence of an active nu- 
cleus ( ,Gonzalez-Delgado et al.|1998,|Colina et al.|2002[ ). [Cid 
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[Fernandes et al. (20011, for a representative sample of 35 
Seyfert 2 galaxies, find that about 40 % of them show un- 
ambiguous evidence of circumnuclear star formation within 
300 pc of the nucleus and that these starforming regions con- 
tribute about 30 to 50 % to the H/3 total emission of the 
central zone. 

These CNSFRs, with sizes going from a few tens to a 
few hundreds of pc (e.g. [Di'az &: Perez-Montero||2000[ ) seem 
to be made of several Hll regions ionised by luminous com- 
pact stellar clusters whose sizes, as measured from high spa- 
tial resolution HST images, are seen to be of only a few pc. 
Their masses, as derived with the use of population syn- 
thesis models in circumnuclear regions of different galaxies, 
suggest that these clusters are gravitationally bound and 
that they might evolve into globular cluster configurations 
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i Maoz et al.||1996 l. The luminosities of CNSFRs are rather 
large with absolute visual magnitudes (M„) between -12 and 
-17 and Ha luminosities between 2 x 10^* and 7 x 10'*'' erg 
s^^ . These values are comparable to those shown by 30 Dor, 
the largest Hll region in the LMC, and overlap with those 
shown by Hii galaxi es ( [Melnick et al.|1988|[Dfaz et al.|2000 



Hoyos fc Dfaz|2006 and references therein) 

Although these Hll regions are very luminous not much 
is known about their kinematics or dynamics for both the 
ionised gas and the stars. It could be said that the worst 
known parameter of these ionising clusters is their mass. 
There are different methods to estimate the mass of a stellar 
cluster. Classically one assumes that the system is virialized 
and determines the total mass inside a radius by applying 
the virial theorem to the observed velocity dispersion of the 
stars (ct*). The stellar velocity dispersion is however hard to 
measure in young stellar clusters (a few Myr old) due to the 
shortage of prominent stellar absorption lines. The optical 
continuum between 3500 and 7000 A shows very few lines 
since the light at these wavelengths is dominated by OB 
stars which have weak absorption lines at the same wave- 
lengths of the nebular emission lines (Balmer If and Hel 
lines). As po inted out by several authors (e.g. |Ho fc Filip 



penko|l996a l, at longer wavelengths (8500 A) the contami- 



nation due to nebular lines is much smaller and since red su- 
pergiant stars, if present, dominate the near-IR light where 
the Can AA8498, 8542, 8662 A triplet hues (CaT) are found, 
these should be easily observable allowing the determination 
of a, ( [Terlevich et al.|1990[|P"rada et al.|19"94l ). We have pre- 



viously detected the CaT lines in CNSFRs but at a spectral 
resolution that was below that required to meeisure accu- 
rately their velocity dispersions (e.g. Terlevich et al.|["l990 1 



Added interest in the study of CNSFRs comes from the 
fact that they are in general of high metal abundance (Diaz 
et al. 2006), therefore they provide clues for the understand- 
ing of star formation phenomena at large metallicities, and, 
being close to the nuclear regions, for the determination of 
metallicity gradients in spiral galaxies. 

is an SBb(r)II spiral galaxy ( |Sandage| 
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by |Sersic fc Pastoriza| ( [T967| . | Alloin fc Nieto| |T9§2l ) studied 
the star formation activity in the nuclear region and along a 
ring of about 20" in diameter using broad and narrow band 
images, concluding that NGC3351 harbours high-mass star 
formation in these zones. In fact it can be considered a nu- 
clear starburst galaxy since the star formation rate per unit 
area in the nuclear region, as estimated from the Hq emis- 
sion, compared to that in the disk is significantly increased 
( Devereux et al.|1992 \ . More recently, near infrared photom- 
etry in the J and K bands has been presented by |Elmegreen| 

et al.| (fl997|l who derive a circumnuclear star formation rate 
1 



of 0.38MQyr-\ From CO emission observations, Planesas 
et al. (1997) derive a mass of molecular gas of 3.5 x 10* M©. 

HST UV images show that the present star formation in 
NGC3351 is distributed along a nuclear ring with a major 
axis of 15" with starforming regions arranged in complexes 
of diameters between 1.4 and 2.0" which are made up of sev- 
eral high surface brightness knots a few pc in size embedded 
in a more diffuse component ( Colina et al.|[l997 |. No signs 
of activity have been observed in the nucleus of NGC 3351. 

There are several values in the literature for the dis- 



tance to NGC 3351. Here we adopt the distance derived by 
Graham et al. (1997) from Cepheid variable stars as part 
of the Hubble Space Telescope (HST) Key Project on the 
Extragalactic Distance Scale, which is 10.05 Mpc giving a 
linear scale of ~49pc per arcsec. 

In this paper we present high-resolution far red spectra 
(~0.39 Apx~^ ~ 13.66kms~^ px~^ at central wavelength, 
Ac=8563A) and stellar velocity dispersion measurements 
(cr,) along the line of sight, for the five circumnuclear 
starforming regions and the nucleus of the spiral galaxy 
NGC 3351. We have also measured the ionised gas ve- 
locity dispersions (cTg) from high-resolution blue spec- 
tra (~0.2lApx-^~ 12.63 kms'Vx"^ at = 4989 A) using 
Balmer H/9 and [Olll] emission lines. The comparison be- 
tween cr, and (Tg might throw some light on the yet unsolved 
issue about the validity of the gravitational hypothesis for 
the origin of the supersonic motions observed in the ionised 
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gas in Giant Hll regions ( Melnick, Tenorio-Tagle fc Terlevich 



Section 2 presents the details of the observations and 
the data reduction. Section 3 presents the results concern- 
ing the kinematics of gas and stars in each of the observed 
regions, as well as the determination of their masses. Sec- 
tion 4 is devoted to the discussion of these results and finally 
the summary and conclusions of this work are presented in 
Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

High resolution blue and far-red spectra were acquired as 
part of an observing run in 2000. They were obtained si- 
multaneously using the blue and red arms of the Intermedi- 
ate dispersion Spectrograph and Imaging System (ISIS), on 
the 4.2 m William Herschel Telescope (WHT) of the Isaac 
Newton Group (ING) at the Roque de los Muchachos Ob- 
servatory, on the Spanish island of La Palma. We used the 
EEV12 detector in the blue arm with a CCD binning factor 
of 2 both in "x" and "y" . The H2400B grating was used to 
cover the wavelength range 4779-5199 A (A^ = 4989 A). The 
TEK4 CCD was attached to the red arm. In this case, the 
R1200R grating was used, covering a spectral range from 
8363 to 8763 A (A^ = 8563 A). The spatial resolution of the 
observations was 0.2 and 0.36" px~^ for the blue and red con- 
figurations respectively. A slit width of 1 arcsec was used. 
The spectral ranges and grating resolutions (spectral dis- 
persion in Apx~^ and FWHM in A measured on the sky 
lines) attained are given in table [l] containing the journal of 
observations. 

Three different slit positions were chosen to observe 5 
CNSFRs and the nucleus of NGC 3351. In Fig.[l]we show the 
position of these slits at angles of 355°, 45° and 310° (which 
we have labelled SI, S2 and S3 respectively) overlapped on 
an image of the circumnuclear region of NGC 3351 acquired 
by the HST (see description below). The regions observed 
are identified on the Ha map by Planesas, Colina fc Perez- 
Olea (1997) using their same nomenclature. The slit was 
positioned on the regions using the offsets from the nucleus, 
clearly visible during the data acquisition, given by Planesas 
et al. (1997). 
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Table 1. Journal of Observations 



Date 


Spectral range 


Disp. 


FWHM 


Spatial res. 


PA 


Exposure Time 




(A) 


(Apx-i) 


(A) 


("px-l) 


n 


(sec) 


04-02-2000 


4779-5199 


0.21 


0.4 


0.38 


355 


4 X 1200 


04-02-2000 


8363-8763 


0.39 


0.7 


0.36 


355 


4 X 1200 


05-02-2000 


4779-5199 


0.21 


0.4 


0.38 


45 


3 X 1200 


05-02-2000 


8363-8763 


0.39 


0.7 


0.36 


45 


3 X 1200 


05-02-2000 


4779-5199 


0.21 


0.4 


0.38 


310 


3 X 1200 


05-02-2000 


8363-8763 


0.39 


0.7 


0.36 


310 


3 X 1200 




Figure 1. Left: F606W (wide V) image centred on NGC 3351 obtained with the WFPC2 camera (PCI) of the Hubble Space Telescope. 
Right: HST-NICMOS image obtained through the F160W filter. For both images the orientation is north up, east to the left. The location 
and P.A. of the WHT-ISIS slit positions, together with identifications of the CNSFRs extracted, are marked. 



Several bias and sky fiat field frames were talten at tlie 
beginning and the end of eacii niglit in botii arms. In addi- 
tion, two lamp flat field and one calibration lamp exposure 
per each telescope position were performed. Tiie calibration 
lamp used was CuNe-|-CuAr. 

We have also downloaded two astrometrically and pho- 
tometrically calibrated broad-band images of the central 
part of NGC 3351 from the Multimission Archive at Space 
TelescopeQ The images were taken through the F606W 
(wide V) and the F160W (H) filters with the Wide Field 
and Planetary Camera 2 (WFPC2; PCI) and the NICMOS 
Camera 3 both on-board the HST. They are displayed in 
Figure 1, with the slit positions overlaid and the CNSFRs 
and the nucleus labelled. 



2.2 Data reduction 

The data was processed and analysed using routines 
in the usual manner. The procedure includes the removal 
of cosmic rays, bias substraction, division by a normalised 
flat field and wavelength calibration. Wavelength fits were 
performed using 20-25 arc lines in the blue and 10-15 lines in 
the far red by a polynomial of second- to third-order. These 
fits have been done at 50 and 60 locations along the slit in 
the blue and far red, respectively, and they have yielded rms 
residuals between ~0.1 and ~0.2px. 

In the red, we also performed the wavelength calibra- 
tion using sky lines following the work by |Osterbrock et al.| 
( |1996| ). However, although more lines (20-25) were available 
for the calibration, the fits gave higher rms residuals, be- 
tween ~0.3 and ~0.45 px since the low intensities of some of 
the lines did not allow a good Gaussian fit. We have there- 
fore adopted the calibration made using only the arc lines. 

Background subtraction was performed using the spec- 
tra at both ends of the slit. This background includes light 



^ http:/ /archive. stsci.edu/hst/wfpc2 



^ IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc. (AURA) under cooperative agreement with the National Sci- 
ence Foundation (NSF). 



4 G. F. Hdgele et al. 



Table 2. Stellar reference frames. 



Star 


ST-LC 


date 


HD71952 


KOI 


04-02-2000 


HD129972 


G6I 




HD134047 


G6III 




HD144063 


G4III 




HD16400 


G5 III 


05-02-2000 


HD22007 


G5I 




HD22156 


G6III 




HD92588 


KlI 




HD102165 


F7I 




HD115004 


G8III 




HD116365 


K3III 





from the disk and bulge of the galaxy. It was almost impossi- 
ble to neatly substract the background bright emission lines 
in the far red spectra due to the nonuniform and extended 
nebular emission surrounding each cluster and the variation 
over time of the sky emission lines. It is worth noting that 
these spurious features do not affect the CaT absorption 
lines. 

We have not corrected the spectra for atmospheric ex- 
tinction or performed any flux calibration, since our purpose 
was to measure radial velocities and velocity dispersions. 

In addition to the galaxy frames, observations of 11 tem- 
plate velocity stars were made (4 during the first night and 7 
during the second) to provide good stellar reference frames 
in the same system as the galaxy spectra for the kinematic 
analysis in the far red. They are late type giant and super- 
giant stars which have strong CaT features (see |Diaz et al.| 
19891. In table [2] we list the spectral types and luminosity 
classes of the stars and the dates of observation. 

Figure [2] shows the spatial profiles of the H/3 emission 
along each slit position. These H/3 profiles have been gener- 
ated by collapsing eleven pixels of the spectra in the direc- 
tion of the resolution between 4860 and 4862 A, rest frame, 
A 4861 being the center of the H/3 line, and are plotted as 
dashed lines. Continuum profiles were generated by collaps- 
ing 11 resolution pixels centered on 4844 A and are plotted 
as dash-dotted lines. The difference betwee the two is shown 
by solid lines and corresponds to the pure H/3 emission. On 
these profiles we have selected the regions of the frames to be 
extracted in one dimensional spectra corresponding to each 
of the identified CNSFRs. Those regions are marked by hor- 
izontal lines and labelled with their respective names in the 
figure. Spectra in slit positions 1 and 2 are placed along the 
circumnuclear ring and therefore any contribution from the 
underlying galaxy bulge is difficult to assess but the spatial 
profiles corresponding to slit position 3, that goes through 
the galactic nucleus, can be used to make an estimate. For 
the blue spectra the light from the underlying bulge is al- 
most negligible amounting at most to 5 % at the H/3 line. For 
the red spectra, the contribution is more important. From 
Gaussian fits to the A 8500 A profile at slit position S3 we 
find that it amounts to about 20 % for the lowest surface 
brightness region, R5. Similar values are found from inspec- 
tion of the HST-NICMOS image shown in Figure 1. 

Figs. |3] and [4] show representative spectra of these re- 
gions (R6 and R3, respectively) in the blue and the red spec- 
tral range split into two panels. The blue spectra show the 
Balmer H/3 recombination line and the collisionally excited 
[Qui] lines at AA 4959, 5007 A. Due to the high metallicity 



of the CNSFRs (Diaz et al. 2006), the lines of oxygen are 
seen to be very weak (see Fig. [5| and, in some cases, only 
the strongest A 5007 A is detected (left panel of this figure). 
The red spectra show the stellar Can triplet (CaT) lines in 
absorption at AA 8498, 8542, 8562 A. In some cases, these 
lines are contaminated by Paschen emission (see for example 
the lower panel of Fig. [4| which occur at wavelengths very 
close to those of the CaT lines. In addition to these emission 
lines, it is possible to observe other emission features, such as 
Oi A 8446, [Chi] A 8579, Pa 14 and [Feii] A 8617. The observed 
red spectrum of R3 is plotted in Fig. |4] with a dashed line. A 
single Gaussian fit to the emission lines was performed and 
the lines were subsequently subtracted after checking that 
the theoretically expected ratio between the Paschen lines 
was satisfied. The solid line shows the subtracted spectrum. 



3 RESULTS 

3.1 Kinematics of stars and ionised gas 

Stars 

Stellar radial velocities and velocity dispersions were ob- 
tained from the absorption CaT lines using the cross- 
correlation technique described in detail by |Tonry fc Dayis] 
( 19791. This method requires the comparison with a stellar 



template (which can be synthetic or observed) that repre- 
sents the stellar population that best reproduces the con- 
spicuous feature used to perform these measurements. An 
example of the red spectrum of a template star (HD 116365) 
is shown in Fig. |6] with the prominent features of CaT in- 
dicated. The line-of-sight velocity dispersions is calculated 
from the width of the primary peak of the cross-correlation 
function (CCF) after deconvolution of the instrumental pro- 
file. A filtering of high frequencies of the Fourier transform 
spectrum is usually included in this procedure to avoid noise 
contamination and a low frequency filtering is usually made 
to eliminate the residual continuum. 

Minor changes/improvements with respect to the cross- 
correlation technique originally proposed by |Tonry fc Davis] 
(11979 1 were introduced as described below. In order to ap- 



ply this method, we have used XCSAO, an external package 
of IRAF within the RVSAO, which implements the cross- 
correlation method of Tonry and Davies and is a direct de- 
scendant of the system built by them (Kurtz fc Mink|1998 1 . 
We used late type giant and supergiant stars that have 
strong CaT absorption lines (Fig. [6| as stellar velocity tem- 
plates. We normalized the stellar spectra dividing by a fit- 
ted continuum and convolved each stellar spectrum template 
with a set of Gaussian functions of different a simulating a 
wide range in velocity dispersions from 10 to 100 km/s with a 
bin size of 5 km/s. The obtained spectra are cross-correlated 
with the original template obtaining a relation between the 
width of the main peak of the cross-correlation and cr of 
the input Gaussian. This relation constitutes a correction 
curve for each template that is applied as described below 
to obtain the stellar velocity dispersion for each CNSFR as 
described in |iSfelson fc Whittle] ( |1995| . This procedure will 
allow us to correct for the known possible mismatches be- 
tween template stars and the region composite spectrum. In 
Figure 7 we show an example of these correction curves, in 
particular for HD 144063, together with a linear fit to it. 
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Figure 2. Spatial profiles of H/3 for each slit corresponding to line+continuum (dashed line), continuum (dashed-dottcd line) and the 
difference between them (solid line), representing the pure emission from H/3. Pixel number increases to the North. Horizontal small lines 
show the location of the CNSFRs and nuclear apertures. 



We determined the line of sight stellar velocity along 
each slit. Extractions were made every two pixels for slit 
position SI. For slit positions S2 and S3 they were made 
every three pixels, with one pixel overlap between consec- 
utive extractions. In this way the signal to noise ratio and 
the spatial resolution were optimised. 

The stellar velocity dispersion was estimated at the po- 
sition of each CNSFR and the nucleus using apertures of 5 
pixels in all cases, which correspond to 1.0" x 1.8". A set of 
eleven templates of different spectral types and luminosity 
classes were used following [Nelson &: Whittie] (|1995| ) with 
the variation introduced by Palacios et al. (19971 of using 



the individual stellar templates instead of an average. To 
measure the velocity dispersion in a galactic spectrum we 
have convolved it with each stellar template, correcting the 
width of the main peak of the CCF (see Fig. |8| with the cor- 
responding correction curve. Although the linear fits to the 
curves are very good approximations (see Fig. [7| , we used 



a linear interpolation between the two nearest values to es- 
timate the corrected width. The a of the stars (cr*) is the 
average of the a values found for each stellar template and 
its error is given by the dispersion of the individual values 
of a and the rms of the residuals of the wavelength fit. This 



procedure allows a more accurate estimate of o", ( Palacios 
et al. 1997). The radial velocities are the average of the ra- 
dial velocities determined directly from the position of the 
main peak of the cross-correlation of each galaxy spectrum 
with each template in the rest frame. 

The stellar velocity dispersions are listed in column 3 
of table [3] along with their corresponding errors. 



Ionised gas 

The wavelength and the width of the H/3 and [Olll] A 5007 A 
emission lines were measured to determine both the radial 
velocities and the velocity dispersions of the ionised gas. As 
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Figure 3. Blue (upper panel) and red (lower panel) rest frame normalised spectra of R6. Notice the absence of conspicuous emission 
lines in the red spectral range for this region. 
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Figure 7. Velocity dispersions correction curve for HD144063 
(circles). The solid line is a linear fit to the curve. 



in the case of the stars, we determined the radial velocity 
of the gas in the line of sight along each slit, every two 
pixels for SI and every three pixels, superposing one pixel 
for consecutive extractions, for S2 and S3. 

The velocity dispersion of the gas was estimated at 
the position of each CNSFR and the nucleus using 5 pixel 
apertures, corresponding to 1.0" x 1.9". Following | Jimenez-] 
Benito et al. ( 2000 1 we adjusted three different suitable con- 




Figure 8. Example of the cross correlation function. 



tinua chosen by visual inspection and fitted a Gaussian to 
the whole line. Positions and widths of the emission lines are 
the average of the corresponding measurements and their er- 
rors are calculated as the dispersion of these measurements 
taking into account the rms of the residuals of the wave- 
length calibration. Thus, the error is associated with the 
continuum placement. 

The velocity dispersions of the gas are calculated as 



where am and ai are the measured and instrumental dis- 
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Figure 4. Blue (upper panel) and red (lower panel) rest frame normalized spectra of R3. The dashed line in the lower panel shows the 
obtained spectrum; the solid line represents the spectrum after subtracting the emission lines (see text). 
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Figure 5. Enlargement of the blue rest frame normalised spectra of R6 (left) and R3 (right). 



persions respectively, ai was measured directly from the sky 
emission lines and is about 10.5 km s~^ at A 4861 A. 

Unexpectedly, the Gaussian fit just described revealed 
the presence of more than one component in the H/3 lines. 
The optimal fit was found for two different components for 
all the regions. We then used the widths of those components 
to fit the [OUT] lines finding also optimal fits for regions 
R2 and R3 in the four analysed spectra corresponding to 
two different slit positions. For the rest of the regions the 
two component fit did not show any improvemt over the 
single component one, probably due to the low signal-to- 



noise ratio. An example of the fit for H/3 and [Olll] A 5007 A 
can be seen in Fig. [9] The radial velocities found using this 
method are the same, within the errors, as those found by 
fitting a single component. 

Gas velocity dispersions for H/3 and [Olll] A 5007 A lines 
are given in table [3] We have listed the velocity dispersions 
derived using the two different methods, the single line and 
the double line Gaussian fit . For the first method we have 
labelled the results as 1 component (columns 4 and 5). For 
the second one, labels 2 components - narrow (columns 6 
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Figure 6. Red rest frame normalised spectrum of HD116365. 
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Figure 9. Sections of the normalised spectrum of R2. The left panel shows from 4865 to 4880 A, containing H/3 and the right panel, 
from 5003 to 5035 A, containing the [Oni] A 5007 A emission line. For both we have superposed the fits from the ngaussfit task in IRAF; 
the dashed-dotted line is the broad component, the dotted line is the narrow component and the dashed line is the sum of both. 



and 7) and 2 components - broad (columns 8 and 9) for the 
narrow and the broad component respectively are used. 



3.2 Sizes 

Two parameters are needed in order to determine the mass 
of a virialized stellar system, namely its velocity dispersion 
and its size. Previously estimated radii (R) for the regions, 
such as those given by Planesas et al. (1997) from Ha images 
or 



Colina et al. (19971 from ultraviolet STIS-HST images, 



were defined to include the total integrated emission flux 
of the regions, and therefore they are not appropriate to 
calculate their dynamical masses, as what is needed for this 
is a measurement of the size of the mass distribution, i.e. 
the star-cluster size. This is traditionally done measuring 
the effective radius using images obtained in bands where 
the light is dominated by the stellar contribution and where 
the contamination by gaseous emission is either small or can 
be estimated and corrected. 



Following Ho & Filippenko 1 1996a I we have defined the 
radius (R) of a given structure as in |Meurer et ah] ( |1995[ ) as- 
suming that they have an intrinsically circularly symmetric 
Gaussian profile. As we can see in Fig. [lO] this is a rather 
good approximation. We have determined the full width at 
half maximum, W50, by superposing the contours in the 



F606W WFPC2-HST image and determining the half light 
strength between the peak of the intensity for each region 
and the galaxy background (red contours in the electronic 
version of the paper). Then, for this model, we take R as the 
half light radius, R = 0.5xW5o- While |Meurer et al.|([l995l 



measured the sizes from an Hq image, we only have the 
F606W (wide V; containing Ha) image in this spectral range 
and with the appropriate spatial resolution. 



We find, as expected, that several of these regions are 
composed of more than one knot. Only R4 and R5 seem to 
have a unique principal knot, at least at the PCI spatial 
resolution. In the most extreme case, R6, we find at least 18 
knots with a detection level of 10 a above background level. 
All these knots are within the radius of the regions defined 
by Planesas et al. (1997). We have to remark that our search 
for knots has not been exhaustive since that is not the aim 
of this work. The radii of the single knots vary between 1.7 
and 4.9 pc. 



Position, radius and its error, and the peak intensity in 
counts measured from the WFPC2 image are listed in table 
[4] for each knot identified. 
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Figure 10. Enlargements of the F606W image around the CNSFRs of our study witli the contours overlapped. The circles correspond 
to the adopted radius for each region. [See the electronic edition of the Journal for a colour version of this figure where the adopted radii 
are in blue and the contours corresponding to the half light brightness are in red.\ 
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Table 3. Velocity dispersions. 



1 component 



2 components 



broad 



Region 


Slit 




0-<,as(H/3) 


tJgas{[Om]) 


(T9as(H/3) 


(Tgas([Om]) 




Ugas{[Olll\) 


R2 


SI 


50±1 


26±1 


72±7 


17±3 


21±4 


45±3 


74±5 


R2 


S3 


51±6 


29±3 


69±9 


16±2 


23±5 


43±2 


76±8 


R3 


SI 


55±5 


35±1 


67±7 


25±3 


28±4 


59±4 


71±4 


R3 


S2 


59±7 


39±5 


70±7 


24±3 


24±6 


59±3 


74±9 


R4 


S2 


66±4 


37±4 


76±8 


29±3 




65±4 




R5 


S3 


47±4 


34±2 


56±7 


30 




76 




R6 


S2 


39±6 


29±6 


46±7 


16±3 




46±4 




N 


S3 


67±1 


53±3 


73±6 


41±5 




67±7 




velocity dispersions in 


km s~^ 













Table 4. Positions, radii and peak intensities. 



Region 


position 


R 


I 




J2000.0 


^.72000. 


(pc) 


(counts) 


R2 


10'"43'°57!12 


+11°42'08'.'27 


4.9±0.2 


2531 


R2a 


10'"43™5rrl2 


+11°42'08'.'53 


2.5±0.2 


956 


R3 


10'i43™57;16 


+11°41'59'.'02 


3.9±0.2 


2598 


R3a 


10'^43""57!16 


+11°41'59'.'24 


2.2±0.1 


1852 


R3b 


10'^43'"57?18 


+11°41'58'.'93 


2.0±0.3 


1032 


R3c 


10'"43™57!21 


+11°41'59'.'41 


2.5±0.3 


546 


R3d 


10'i43™57;21 


+11°41'59'.'33 


2.7±0.3 


489 


R3e 


10'^43""57!20 


+11°41'59'.'71 


2.9±0.2 


383 


R3f 


10'^43™57?20 


+11°41'59'.'26 


3.7±0.4 


330 


R4 


10'"43™57?32 


+11°42'00'.'92 


2.9±0.1 


2581 


R5 


10'"43™57!55 


+11°42'02'.'50 


3.2±0.1 


973 


R6 


10'"43™57!60 


+11°42'05'.'21 


2.7±0.1 


2399 


R6a 


10'i43™57?62 


+11°42'06'.'81 


2.0±0.1 


794 


R6b 


10'^43""57!62 


+11°42'06'.'00 


2.0±0.2 


714 


R6c 


10'^43'"57?63 


+11°42'06'.'38 


2.0±0.1 


551 


R6d 


10'"43™57!62 


+11°42'05'.'86 


1.7±0.1 


404 


R6e 


10'i43'"57:65 


+11°42'06'.'85 


2.4±0.2 


376 


R6f 


10'^43'"57!62 


+11°42'06'.'21 


2.0±0.1 


340 


R6g 


10'^43™57?61 


+11°42'06'.'12 


1.7±0.1 


295 


R6h 


10'^43™57!63 


+11°42'06'.'52 


1.7±0.1 


290 


R6i 


10'^43™57?64 


+11°42'06'.'13 


2.0±0.1 


285 


R6j 


10'^43""57!62 


+11°42'05'.'66 


2.9±0.3 


282 


R6k 


10'^43""57?61 


+11°42'05'.'41 


2.4±0.2 


275 


R61 


10'"43™57?61 


+11°42'06'.'31 


2.4±0.3 


272 


R6m 


10'i43™57:61 


+11°42'05'.'43 


2.2±0.2 


270 


R6n 


10'^43™57!63 


+11°42'05'.'52 


2.2±0.3 


251 


R6o 


10'^43'"57!63 


+11°42'04'.'83 


2.9±0.3 


212 


R6p 


10'"43™57!61 


+11°42'04'.'81 


3.1±0.3 


212 


R6q 


10'"43'"57?56 


+11°42'05'.'48 


2.2±0.3 


211 


N 


10'i43™57?29 


+11°42'05'.'81 


11.3±0.1 


513 



3.3 Masses 

Dynamical masses 

We have used the virial theorem to estimate upper limits 
to the dynamical masses (M*) inside the half light radius 
(R) for each observed knot in the F606W WFPC2-HST im- 
age. In order to do this, we have assumed that the systems 
are spherically symmetric, gravitationally bound and have 
isotropic velocity distribution (cr'^ (total) = 3 a ^). Then, the 



dynamical mass is given by M. = 3 at R/G ( |Ho fc Filip- 
penko|1996a|b| ). 

It must be noted that we have measurements for the size 
of each knot (typically 5pc), but we do not have direct ac- 
cess to the stellar velocity dispersion of each of the clusters. 



since our spectroscopic measurements encompass a wider 
area (1.0" x 1.8" ~ 49 x 88 pc^) that includes the CNSFRs 
to which each group of knots belong. The use of these wider 
size scale velocity dispersion measurements to estimate the 
mass of each knot, lead us to overestimate the mass of the 
individual clusters, and hence of each CNSFR. 

However, as can be seen in the HST-NICMOS image 
(right panel of Fig. [l]), the CNSFRs of NGC 3351 are clearly 
visible in the IR and dominate the hght inside the apertures 
observed. All the regions analysed show up very prominently 
in the near IR and therefore we can assume that the hght 
at the CaT wavelength region is dominated by the stars in 
the clusters. The IR Call triplet is very strong, in fact the 
strongest stellar feature, in very young clusters, i.e. older 



than 4 Myr ( Terlevich et al. 1990 1 . Besides, we detect a mini- 



mum in the velocity dispersion at the position of the clusters, 
indicating that they are kinematically distinct. We cannot 
be sure though that we are actually measuring their veloc- 
ity dispersion and thus prefer to say that our measurements 
of CT* and hence dynamical masses constitute upper limits. 
Although we are well aware of the difficulties, still we are 
confident that these upper limits are valid and important 
for comparison with the gas kinematical measurements. 

The estimated dynamical masses and their correspond- 
ing errors for each knot are listed in table |5] For the regions 
that have been observed in more than one slit position, we 
list the derived values using the two different stellar veloc- 
ity dispersions. The dynamical masses in the lines labelled 
"sum" are the sum of the values derived for each knot of the 
corresponding region, when there are more than one. The 
"adopted" dynamical mass for R2 is directly the sum of the 
dynamical masses of the two knots derived using the velocity 
dispersion estimated from SI since the mass calculated from 
S3 is in agreement with the first one (see table [sj but it has 
a much greater error. This greater error is due to the lower 
signal to noise ratio of the spectrum of R2 extracted from 
S3. In contrast, the derived dynamical masses for R3 are 
not in complete agreement, although they coincide within 
the errors, which are comparable. In this case we have taken 
the average as the "adopted" value. The fractional error in 
the "sum" and the adopted values of the dynamical masses 
of the CNSFRs are listed in column 4 of table (5) 



Ionising star clusters 

We have derived the masses of the ionising star clusters 
(Mion) from the total number of ionising photons using so- 
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Table 5. Dynamical masses. 



Region 


Slit 


Mt 


error (%) 


R2 


SI 


85±5 




R2a 


SI 


44±4 




R2sum 


SI 


129±6 


3.4 


R9 


oo 






R2a 


S3 


44±11 




R2sum 


S3 


131±23 


17.8 


Kz (aaoptca j 




129it6 


3.4 


R3 


SI 


82±14 




R3a 


SI 


46±8 




R3b 


SI 


42±9 




R3c 


SI 


52±11 




R3d 


SI 


57±11 




tXOG 


O i 


-1-1 1 

Dim 1 ± 




R3f 


SI 


78±15 




R3sum 


SI 


417±31 


7.3 


R3 


S2 


94±23 




R3a 


S2 


53±13 




R3b 


S2 


48±14 




R3c 


S2 


60±16 




R3d 


S2 


65±17 




R3e 


S2 


70±17 




R3f 


S2 


89±5 




R3sum 


S2 


477±47 


9.9 


R3 (adopted) 




447±56 


12.6 


R4 


S2 


87±12 


13.9 


R5 


S3 


49±8 


16.4 


R6 


S2 


29±9 




R6a 


S2 


21±8 




R6b 


S2 


21±8 




R6c 


S2 


21±8 




R6d 


S2 


18±7 




R6e 


S2 


26±10 




R6f 


S2 


21±8 




R6g 


S2 


18±7 




R6h 


S2 


18it7 




R6i 


S2 


21±8 




R6j 


S2 


31±12 




R6k 


S2 


26±10 




R61 


S2 


26±10 




R6m 


S2 


24±9 




R6n 


S2 


24±9 




R6o 


S2 


31±12 




R6p 


S2 


33±13 




R6q 


S2 


24±9 




R6sum 


S2 


434±39 


9.0 


N 


S3 


350±11 


3.1 


masses in 10^ M 










lar metallicity single burst mode ls by [Garcia- Vargas et al.| 
( 1995 1, assuming a Salpeter IMF ( Salpeter|1955 l with lower 
and upper mass limits of 0.8 and 120 Mq which provide the 
number of ionising photons per unit mass, (Q{Ho) / Mion) ■ 
This number decreases with the age of the region. We have 
used the following relation between Q{Ho)/Mion and the 
equivalent width of H/3, EW(H/?), derived from the models, 
in order to take into account the evolutionary state of the 



region I Dfaz||1998j) : 



We have taken the values of EW(H/3) from Perez-Olea 
( 1996 1 (see table|6| which are not corrected by the contribu- 
tion to the continuum by the underlying stellar population. 
This correction would increase the values of EW(H/3) thus 
decreasing the calculated Mion- The total number of ion- 
ising photons has been derived from the Ha luminosities 
( Leitherer fc Heckman|1995 1: 



Q{Ho) = 7.35 X 10"L(i/Q) 

We have taken the total observed Ha luminosities from 
Planesas et al. (1997) correcting them for the different as- 
sumed distance. These authors estimated a diameter of 2.4" 
for each whole region and the nucleus, except for R4 for 
which they used 2.2". We have corrected the Ha luminosi- 
ties for internal extinction using the colour excess [E(B-V)] 
estimated by Pcrcz-01e a| (|1996[) and assuming the galactic 
extinction law of Mi ller fc Mathews] ( [1972) with i?„ = 3.2. 

Our derived values of Q(Ho) are lower limits since we 
have not taken into account either the absorption of photons 
by dust or any photon escape from the Hll regions. 

The final expression for the derivation of Mion is: 



Mion = 



7.35 X 10"L(//a) 

IQ44.8 + 0.86 loglEW(H0)] 



Ionised gas 

The amount of ionised gas (Mhii) associated to each star- 
forming region complex has been obtained from our derived 
Hq lu minosities using the relation given by |Macchetto et al.| 
( 1990 1 for an electron temperature of 10* K 



Mh 



= 3.32 X 10"^^ L(ifa) iV; 



log (q (Ho) /Mion) 



44.8 -f- 



0.86 log (^EW{H 13)) 



where Ne is the electron density. The electron density for 
each region (obtained from the [Sll] AA 6717 / 6731 A line ra- 
tio) has been taken from Di'az et al. (2006) for the CNSFRs 
and Perez-Olea (19961 for the nucleus (see table |6|. 

In table |6| we have listed for each region and the nu- 
cleus, the Hq luminosities [L(Ha)] (corrected and uncor- 
rected for reddening), the colour excess, E(B-V), the log- 
arithmic extinction at Ha, C(Ha), the number of ionising 
photons, Q(Ho), the equivalent width of H/3, EW(H/3), the 
electron density, Ne, and the masses of the ionising stellar 
cluster, Mion, and of the ionised hydrogen, Mhii- 



3.4 Emission line ratios 

We have used two different ways to integrate the intensity 
of a given line: (1) in the cases of a single Gaussian fit the 
emission line intensities were measured using the SPLOT 
task of IRAF. The positions of the local continua are placed 
by eye. For the H/3 emission lines a conspicuous underly- 
ing stellar population is easily appreciable by the presence 
of absorption features that depress the lines (see discussion 
in |Dfaz 1988[ ). Examples of this effect can be appreciated 
in Fig. |5| We have then defined a pseudo-continuum at the 
base of the line to measure the line intensities and minimise 
the errors introduced by the underlying population (for de- 
tails see Hagele et al.||20d6 l. (2) in the cases of a fit by two 
Gaussians the individual intensities of the narrow and broad 
components are estimated from the fitting parameters using 
the expression I = 1.0645 A x FWHM (= y/2n A x a), where 
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Table 6. Physical parameters. 



Region Lof,s(Ha)t 


E{B-V)t 


C(Ha) 


L(Ha) 


Q(Ho) 


EW(H/3)t 




Ni 


Mhii 


M,„„/M,(%) 


R2 


19.3 


0.17 


0.02 


20.5 


15.0 


9.5 


7.2 


440 


0.15 


5.6 


R3 


25.0 


0.46 


0.65 


111.0 


81.5 


16.5 


24.2 


430 


0.86 


5.4 


R4 


12.7 


0.27 


0.24 


22.0 


16.2 


13.0 


5.9 


310 


0.24 


6.8 


R5 


5.9 


0.25 


0.20 


9.2 


6.8 


5.1 


5.5 


360 


0.09 


11.3 


R6 


7.5 


0.0 


-0.34 


3.4 


2.5 


2.3 


4.1 


360 


0.03 


1.0 


N 


3.3 


0.07 


-0.19 


2.1 


1.6 


1.8 


3.1 


650* 


0.01 


0.9 


luminosities 


in lO-^** erj 


;s^^, masses 


in 10^ 


M0, ionisin 


g photons 


in 10™ phs" 


^ and densities 


in cm" 


-3 



tfrom Planesas et al. (1997) corrected for the different adopted distances 
tfrom Perez-Olea (1996) 
§ from Diaz et al. (2006) 



I is the Gaussian intensity, A is the amplitude of the Gaus- 
sian and FWHM is the full width at half maximum {a is the 
dispersion of the Gaussian). A pseudo-continuum for the H/3 
emission line was also defined in these cases. 



Following Gonzalez-Delgado et al. (19941, Castellanos, 



Di'az & Terlevich (2002) and Perez-Montero & Di'az (20031, 



the statistical errors associated with the observed emis- 
sion fluxes have been calculated using the expression 
ai=acN^^^[l + EW/(NA)]^/^, where ai is the error in the 
observed line flux, CTc represents the standard deviation in 
a box near the measured emission line and stands for the 
error in the continuum placement, N is the number of pixels 
used in the measurement of the line intensity, EW is the 
line equivalent width, and A is the wavelength dispersion in 
angstroms per pixel. For the If/3 emission line we have dou- 
bled the derived error, cr; , in order to take into account the 
uncertainties introduced by the presence of the underlying 
stellar population ( Hagele et al.|[2006 K 

The logarithmic ratio between the emission line inten- 
sities of [Olll] A 5007 A and H/3, and their corresponding er- 
rors, are presented in table [7] We have also listed in this 
table the logarithmic ratio between the emission line fluxes 
of [Nil] A 6584 A and Ha together with their corresponding 
errors from Di'az et al. (20061. 



4 DISCUSSION 

The upper panel of Fig. [TT] shows the relation between the 
velocity dispersions of gas and stars. Circles correspond to 
gas velocity dispersions measured from the H/3 emission line 
using a single Gaussian fit. Squares and triangles correspond 
to measurements performed using two-component Gaussian 
fits, squares for the broad component and triangles for the 
narrow one. The straight line shows the one-to-one relation. 
As a general result, the H/3 velocity dispersions derived by a 
single Gaussian fit are lower than the stellar ones by about 
25 kms~^. This is also the case for the Paschen lines in the 
two regions where they could be measured. On the other 
hand, a good agreement is found between the stellar veloc- 
ity dispersions and those of the broad component of H/3. The 
deviant point, marked with arrows, corresponds to the re- 
gion with the lowest signal-to-noise ratio (R5) for which the 
fits do not provide accurate results. The narrow component 
shows velocity dispersions even lower than those obtained 
by single Gaussian fits. The ratio between the fluxes in the 
narrow and broad components is between 0.7 and 0.95 (ex- 



cept for the case of R5 for which no meaningful result is 
found). 

The lower panel of Fig. [TT] shows the relation between 
the stellar velocity dispersions and those of the [Olll] emis- 
sion line measured by both single and two-component Gaus- 
sian fits, which due to the weakness of the line (see figures 
3 to 5), could be done in only four cases, corresponding to 
the two slit positions on regions R2 and R3. In this case, 
the broad component seems to dominate the width of the 
emission line which again agrees with the stellar one. The 
width of the narrow component, in the cases in which the 
two-component fit was possible, is comparable to that of the 
narrow component of the H/3 line. 

The two gaseous components have been plotted sepa- 



rately in Fig. 12 which shows their location in the [Olll]/H/3 



vs [N il] /Ha diagram (BPT; [Baldwin, Phillips fc Terlevich 
[1981[ ) together with a sample of emission line galaxies (in- 
cluding Hll-like objects) from the S DSS-DR3 (Lopez. 2005[) 
and Hll regions from the sample of |Perez-Montero fc Diaz[ 
( |2005[ ). The two systems are clearly segregated in this di- 
agnostic diagram with the narrow component showing the 
lowest excitation of the two and occupying the same posi- 
tion in the diagram as the starburst systems in the SDSS 
dataset with the lowest excitation found. 

Fig.[T3[shows the radial velocities along the slit for each 
slit position as derived from the ionised gas emission lines 
and stellar absorptions. The turnover points of the rota- 
tional curves seem to be located at the same position than 
the starforming ring, specially for the S3 slit position that 
crosses the centre of the galaxy, as found in other galaxies 
(see Telesco fc Decher|1988 Diaz et al.|1999 and references 
therein). For the systemic velocity of NGC 3351, the derived 
values are consistent with those previously obtained by |Ru-| 
bin et al. (19751 and Planesas et al. (1997), and with the 



velocity distribution expected for this type of galaxies ( Bin- 
ney fc Tremaine|1987[ ). 



The rotation velocities derived for both stars and gas 
are in reasonable agreement, although in some cases the gas 
seems to rotate somewhat faster than the stars. In fact, in 
the lower panel of Fig 13, which corresponds to the slit po- 
sition passing through the nucleus, it is interesting to note 
that at the maximum and minimum of the velocity curve, 
which correspond approximately to the positions of regions 
R5 and R2, the H/3 emission line shows velocities lower and 
higher than the stars by about 25 and 30kms~^ respec- 
tively. This could be interpreted as motions of the ionised 
hydrogen deviating from rotation and consistent with a ra- 
dial infall to the central regions of the galaxy. A similar 
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Table 7. Line ratios. 







1 component 


2 components 










narrow 


broad 




Region 


Slit 


log([Olll]5007/H/3) 


log([Olll]5007/H/3) 


log([Oni]5007/H/3) 


log([Nn]6584/H/3)t 


R2 


SI 


-1.07±0.06 


-1.66it0.08 


-0.93it0.12 


-0.43it0.01 


R2 


S3 


-I.01it0.06 


-1.55it0.08 


-0.96it0.13 




R3 


SI 


-I.10it0.06 


-1.57it0.07 


-0.93it0.10 


-0.42it0.01 


R3 


S2 


-1.00it0.06 


-1.52it0.09 


-0.89it0.10 




R4 


S2 


-1.03it0.07 






-0.49itD.01 


R5 


S3 


-0.85it0.12 






-0.37it0.03 


R6 


S2 


-1.09it0.11 






-0.52it0.02 


N 


S3 


-0.28it0.05 









Tfrom Di'az et al. (2006) 



result was found by Rubin et al. (1975) from high disper- 
sion observations of the Ha line in the central region of this 
galaxy. Their preferred model for the fitting of the kinemat- 
ical data consists of gas which is rotating and contracting 
with Vrot = 126 ± 16 km s^^ and Ncont = 34 ± 11 kms^^. 

It is clearly of major interest to find out how widespread 
is the presence of two distinct components in the emission 
lines in ionised regions and what is its influence on the ob- 
served line ratios for several reasons. Firstly, a change in po- 
sition in the diagnostic diagrams would certainly affect the 
classification of the activity in the central regions of the con- 
cerned galaxies. Secondly, it will affect the inferences about 
the nature of the source of ionisation in the two components. 
And thirdly, it could have an influence on the gas abundance 
determinations given that the ratio of the narrow to the 
broad components of H/3 is about 1. Clearly it is not possi- 
ble to use global line ratios to estimate gaseous abundances 
if the permitted and forbidden line fluxes are partially orig- 
inated in different kinematical systems. To disentangle the 
origin of these two components it will be necessary to map 
these regions with high spectral and spatial resolution and 
much better S/N in particular for the O^"*" lines. 3D spec- 
troscopy with Integral Field Units (IFUs) would be the ideal 
tool to approach this issue. 

Our values of the dynamical masses have been derived 
from stellar velocity dispersion measurements as mapped 
by the CaT absorption lines and the sizes measured on an 
HST image. They are in the range between 4.9 x l(f and 
4.3 X 10^ M0 for the CNSFRs and is 3.5 x lO'^ Mq for the 
nuclear region inside the inner 11.3 pc (see table [5|. The 
fractional errors of the dynamical masses of the CNSFRs are 
between ~3.4 and ~ 17.8 %, and is ~3.1 % for the nucleus. 
Masses derived from the H/3 velocity dispersion under the 
assumption of a single component for the gas would have 
been underestimated by factors between approximately 2 
and 4. We have found that only two of the regions, R4 and 
R5, seem to possess just one knot showing up in the con- 
tinuum image and coincident with the Ha emission. The 
sizes of these knots are 2.9 and 3.2 pc respectively. The rest 
of the regions are made up of several knots which presum- 
ably correspond to individual star clusters. Their sizes are 
between 4.9 and 1.7 pc, at the limit of the resolution. For 
comparison, the size of cluster A in NGC 1569 is 1.9 pc, as 



ranges established bylMeurer et al. (19951 in the definition 



given by Meurer et al. ( 1995 1. The sizes and absolute visual 



of a super star cluster (SSC). 

Other estimates of the masses of these regions have 
been obtained from the fitting of broad band colours or 
spectra with the use of stellar population synthesis mod- 
els. Elmegreen et al. (1997) from near IR photometry in 
the J and K bands and models by Leitherer & Heckman 
(1995) for instantaneous star formation and solar and twice 
solar metallicity, derive masses of the CNSFRs from 1 to 
IOxIO^Mq. Colina et al. (1997), from UV (lUE) spec- 
tra an d instantaneous burst models by |Mas-Hesse fc Kunth| 



( 1991 1 derive a value of 3 x lO'^ Mq for the whole SF ring. In 



both cases, in fact, these observables trace the young massive 
population which constitutes only part of the total mass. 

The masses derived here for the circumnuclear starform- 
ing individual clusters of NGC 3351 are between 1.8 and 
8.7 X 10^ Mq. These values are between 5.5 and 26 times 
the mass deriv ed for the SSC A in NGC 1569 by |Ho fc Filip^ 
penko ( 1996a I from stellar velocity dispersion measurements 
using red (~6000 A) spectra, which is (3.3±0.5) x lO'^ Mq 
and larger than masses derived kinematically for SSC in ir- 
regular galaxies fM cCrady, Gilbert fc Graham|2003 Larsen, 
Brodie & Hunter 20041. They are also larger than those de- 



rived by [Benedict et ai.| ( |2002[ ) for the individual circumnu- 
clear clusters of NGC 4314 form HST imaging data following 



the procedure of Elmegreen & Salzer (19991, which are in 
the range 0.2 x 10* < M< 1.6x 10* Mq. However, also the 
H/3 luminosities of the NGC 4314 clusters are lower than 
those of NGC 3351 by a factor of about 40. 

The masses of the ionising stellar clusters of the CNS- 
FRs, as derived from their Ha luminosities, are between 
4.1 X 10^ and 2.4 x 10® Mq for the starforming regions, and 
is 3.1 X 10^ Mq for the nucleus (see table [6|. These values 
are comparable to that derived by Gonzalez-Delgado et al. 
(1994) for the circumnuclear region A in NGC 7714 (5.1 x lO'' 
Mq). In column 11 of table[6]we show a comparison (in per- 
centage) between ionising stellar masses of the circumnu- 
clear regions and their dynamical masses. These values vary 
approximately in the range 1-11% for the CNSFRs, and is 
0.9 % for the nucleus. Since the CaT absorption features are 



magnitudes estimated for each individual star cluster (e.g. 
M„ =-12.55 for R4 derived from the HST image) are in the 



dominated by young stars (see discussion above, Sec. 3.31, 
and the Mion/M* fraction is still remarkably small in the 
case of the CNSFRs composed of single knots (R4 and R5, 
for which the dynamical mass is most robustly estimated) we 
can assume that our upper limits to the dynamical masses, 
in spite of the limitations of the method used to derive them. 
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Figure 12. BPT diagram [Olll]/H/3 vs [Nll]/Ha. Squares corre- 
spond to tlie ratio of the emission intensities of [Olll] and H/3 
estimated using a single Gaussian fit, diamonds to the narrow 
components of the two Gaussian fits and circles to the broad 
components. Dotted and dashed curves are the boundary between 
Active Galactic Nuclei (AGNs) and Hll galaxies defined by Kew- 
ley et al. (200f) and Kauffmann et al. (2003) respectively. The 
solid horizontal and vertical lines represent the division between 
Seyfert galaxies and LINERs according to Ho et al. {f997). Dots 
correspond to a subsample of emission line objects, including Hll 
galaxies, from SDSS-DR3, from Lopez (2005). Triangles corre- 
spond to Hll regions from the sample of Perez-Montero & Diaz 
(2005). They have been split into low metallicity (upside down 
triangles) and high metallicity regions (upward triangles) accord- 
ing to the criterion by Diaz & Perez-Montero (2000) based on 
oxygen and sulphur abundance parameters. 



kinematics follows that of the stars in the regions, all derived 
quantities would be smaller by a factor of two. 



Figure 11. Upper panel: relation between velocity dispersions of 
the gas (derived from H/9) and stars (CaT) for the CNSFRs and 
the nucleus. Symbols are as follows: single Gaussian fit, circles; 
two Gaussian fit, broad component: squares; narrow component, 
triangles. Lower panel, as the upper panel for the [Gin] line. Up- 
ward triangles correspond to the estimates using a single Gaussian 
fit, diamonds represent the broad components of the two Gaussian 
fit and left triangles, the narrow components. 



are rather tight. Then, our results concerning the Mio„/M, 
fraction are robust. 

Finally, the masses of the ionised gas vary between 
3 X 10^ and 8.6 x lO" Mq for the CNSFRs, and is 1 x 10^ Mq 
for the nucleus (see table [6|, also comparable to that derived 
by Gonzalez-Delgado et al. (1994) for the circumnuclear re- 
gion A in NGC7714 (3x10^ Mq). They make up a smaU 
fraction of the total mass of the regions. 

Both the masses of the ionising stellar clusters and the 
ionised gas, have been derived from the Ha luminosity of 
the CNSFRs assumed to consist of one single component. 
However, if we consider only the broad component whose 



5 SUMMARY AND CONCLUSIONS 

We have measured gas and stellar velocity dispersions in 5 
CNSFRs and the nucleus of the barred spiral NGC 3351. The 
stellar dispersions have been measured from high resolution 
spectra of the CaT hues at AA 8494, 8542, 8662 A, while the 
gas velocity dispersions have been measured by Gaussian fits 
to the H/3 A 4861 A line on high dispersion spectra. 

Stellar velocity dispersions are between 39 and 
67kms~^, about 20kms~^ larger than those measured for 
the gas. However, the best Gaussian fits involved two dif- 
ferent components for the gas: a "broad component" with 
a velocity dispersion similar to that measured for the stars, 
and a "narrow component" with a dispersion lower than the 
stellar one by about 30kms~^. 

When plotted in a [Olll]/H/3 vs [Nii]/Hq: diagram, the 
two systems are clearly segregated with the narrow compo- 
nent having the lowest excitation and being among the low- 
est excitation line ratios detected within the SDSS dataset 
of starburst systems. 
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Figure 13. Radial velocities along the slit vs pixel number for 
each slit position (upper panel: SI; middle panel; S2; lower panel: 
S3) as derived from the gas emission lines (circles: H/3; triangles: 
[Olll]) and the stellar absorption ones (squares). The individual 
CNSFRs and the nucleus, "N", or the closest position to it, "d", 
are marked in the plots. The dashed-dotted line is the systemic 
velocity of NGC 3351 derived by Planesas et al. (1997). The dis- 
tance in arcsec from the nucleus is displayed in the upper x-axis 
of each panel. 



Values for the dynamical masses of the CNSFRs have 
been derived from stellar velocity dispersions and are in the 
range between 4.9 x lO'' and 4.3 x 10^ Mq for the CNSFRs 
and is 3.5 x 10^ M© for the nuclear region inside the inner 
11.3 pc. Masses derived from the H/? velocity dispersion un- 
der the assumption of a single component for the gas would 
have been underestimated by factors between approximately 
2 to 4. 

The derived masses for the individual clusters of 
NGC 3351 are between 1.8 and 8.7 x 10*^ Mq. These values 
are between 5.5 and 26 times the mass derived for the SSC 



A in NGC 1569 by Ho fc Filippenko (1996a I and larger than 



other kinematically derived SSC masses. 

Masses of the ionising stellar clusters of the CNSFRs 
have been derived from their Ha luminosities under the as- 
sumption that the regions are ionisation bound and without 
taking into account any photon absorption by dust. Their 
values are between 4.1 x 10^ and 2.4 x 10® Mq for the star- 
forming regions, and is 3.1 x 10^ M© for the nucleus (see 
table [bJi, comparable to that derived by Gonzalez-Delgado 
et al. (1994) for the circumnuclear region A in NGC 7714. 
Therefore, the ratio of the ionising stellar population to the 
total dynamical mass is between 0.01 and 0.11. 

Derived masses for the ionised gas vary between 3 x 10'^ 
and 8.6 x 10^ Mq for the CNSFRs, and is 1 x 10^ Mq for 
the nucleus, also comparable to that derived by Gonzalez- 
Delgado et al. (1994). 

It is interesting to note that, according to our findings, 
the SSC in CNSFRs seem to contain composite stellar pop- 
ulations. Although the youngest one dominates the UV light 
and is responsible for the gas ionisation, it constitutes only 
about 10 % of the total. This can explain the low equiva- 
lent widths of emission lines measured in these regions. This 
may well apply to the case of other SSC and therefore con- 
clusions drawn from fits of SSP (single stellar population) 



models should be taken with caution (e.g. McCrady, Gilbert 
fc Graham]|2003[ [Larsen, Brodie fc Hunter||2004[ ). Also the 
composite nature of the CNSFRs means that star formation 
in the rings is a process that has taken place over time peri- 
ods much longer than those implied by the properties of the 
ionised gas. 

The rotation velocities derived for both stars and gas 
are in reasonable agreement, although in some cases the 
gas shows a velocity slightly different from that of the 
stars. The rotation curve corresponding to the position 
going through the centre of the galaxy shows maximum and 
minimum values at the position of the circumnuclear ring, 
as observed in other galaxies with CNSFRs. The differences 
in velocity between gas and stars can be interpreted as 
motions of the ionised hydrogen deviating from rotation 
and consistent with a radial infall to the central regions of 
the galaxy. Our results are consistent with those found by 
Rubin et al. (1975) and would yield an infall velocity of 
about 25kms~^. 

The existence of more than one velocity component 
in the ionised gas corresponding to kinematically distinct 
systems, deserves further study. Several results derived 
from the observations of the different emission lines could 
be affected, among others: the classification of the activity 
in the central regions of galaxies, the inferences about 
the nature of the source of ionisation, the gas abundance 
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determinations, the number of ionising photons from a 
given region and any quantity derived from them, etc. 
To disentangle the origin of these two components it will 
be necessary to map these regions with high spectral and 
spatial resolution and much bettor S/N in particular for 
the O^"*" lines. 3D spectroscopy with Integral Field Units 
(IFUs) would be the ideal tool to approach this issue. 
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